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Predictive  Flow  Control  to  Minimize  Convective  Time  Delays 
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O0  overcome  the  convective  time  delay  issue  for  active  closed-loop  flow  control  a  model  based  predictive 
control  algorithm  is  analyzed.  From  a  controls  perspective,  forms  of  internal  model  control  or  model  predictive 
control  can  be  used  to  accommodate  and  minimize  the  effect  of  systems  w-ith  pure  time  delays.  Moreover,  the 
Smith  predictor  is  a  commonly  employed  control  technique  to  negate  the  pure  time  delay  in  a  closed-loop  system. 
This  form  of  model  predictive  control  is  applied  to  the  asymmetric  vortex  problem  of  an  axisymmetrie  forebody 
(specifically  a  von  Kerman  ogive  with  fineness  ratio  of  3,5)  at  high  angle  of  attack  (a  —  50*)*  Full  order  Navier- 
Stokes  are  numerically  solved  and  provide  the  plant  process.  Small  port  and  starboard  blowing  patches  are  used  to 
introduce  fluidic  disturbances  at  the  nose  of  the  ogive  to  augment  the  global  flow  stale  and  produce  a  deterministic 
vortex  stale.  As  the  active  flow  control  technique  exploits  ihc  convective  instability,  a  convective  time  delay  exists, 
Linear-timednvariant  and  non-linear  time  invariant  models  are  developed  from  the  open- loop  dynamics.  A  Smith 
predictor  is  employed  within  the  full  order  CFD  simulation.  The  results  of  the  predictive  control  arc  compared 
to  open-loop  and  model- free  closed  loop  behavior.  It  is  shown  that  the  predictive  control  developed  in  this  paper 
while  very  suitable  for  control  of  this  type  of  flow  is  very  sensitive  to  modeling  uncertainties. 


Nomenclature 


A 

Base  area 

V 

Velocity  vector 

c 

Linear  observation  matrix 

P 

Pressure 

Cy 

Estimated  *sidc  force  coefficient 

Re 

Reynolds  number*  UooD/v 

Cp 

Pressure  coefficient 

t 

Time 

D 

Diameter 

V 

Freestream  velocity 

D 

Diameter 

Ul,U2 

Reduced-order  model  inputs 

D 

Base  diameter 

Uoo 

Free-stream  velocity 

/, 

Fineness  ratio,  Lt:Ctnft/D 

X* 

Surface  pressure  sensor  location 

L 

Length  of  the  ogive  model 

L 

Model  length 

a 

Angle  of  attack 

Uaft 

Aft  body  length 

A  CP 

Differential  pressure  coefficient,  port  -  star¬ 

■^WTie 

Nose  cone  length 

board 

M 

Mach  number 

0 

Circumferential  angular  position 

Kp,Tz 

.  C  Tu>  Td  Reduced-order  model  parameters 

P 

Density 

n 

Normal  vector 

T 

Convective  time 

I.  Introduction 

Research  in  dosed- loop  flow  control  is  attempting  to  improve  the  performance  and  robustness  of  control 
strategies  for  highly*  non-linear  complex  flow  fields.  Due  to  natural  instabilities  in  fluid  flows*  closed- loop 
flow  control  has  the  potential  to  improve  performance  of  all  engineering  systems  in  which  a  fluid  flows,  e,g. 
external  flows  around  air  vehicles  or  ground  based  systems  such  as  bridges  and  buildings*  internal  flows  in  pipes 
and  propulsion  systems,  acoustical  emission  and  mitigation*  combustion  and  mixing  problems,  turbulence  and 
transition  management,  as  well  as  alternative  energy  applications  such  as  wave  and  wind  power  devices.  Flow 
control  can  be  categorized  into  two  approaches,  passive  and  active.  Passive  flow  control  makes  use  of  geometry 
modifications  to  alter  the  mean  fluidic  behavior  to  produce  a  desired  effect.  A  classic  example  of  passive  control 
is  the  dimples  on  a  golf  ball  to  transition  the  boundary  layer,  thus  delaying  separation  and  reducing  pressure  drag. 
Another  example  from  airplane  aerodynamics  is  the  boundary'  layer  fences  on  wings  that  ensure  that  longitudinal 
vortices  are  created  and  delay  flow  separation  at  high  angles  of  attack.  Such  passive  flow  control  strategics  arc 
based  on  a  quasi -steady  flow  assumption  and  have  been  extensively  studied  in  the  past. 

In  contrast,  active  flow  control,  which  introduces  a  small  amount  of  energy  to  the  flow  (e,g,  momentum  addi¬ 
tion,  mass  addition,  vorticity  addition),  allows  for  direct*  local,  and  instantaneous  modification  of  the  flow  field. 
The  most  efficient  way  to  employ  active  flow  control  is  to  exploit  natural  instabilities  in  the  flow  field,  which  al¬ 
lows  for  a  minimal  actuation  effort.  A  perturbation  is  introduced  into  the  flow  and  is  amplified  through  the  fluidic 
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instability  as  the  disturbance  convects  downstream  in  the  flow  field:  this  phenomenon  is  commonly  referred  to 
as  a  convective  instability.  Analyzing  published  active  flow  control  efforts,  a  number  of  technical  challenges  can 
be  identified  In  large,  the  temporal  delay  due  to  the  time  necessary  for  a  perturbation  to  convert  and  grow'  to  a 
significant  size  to  effect  the  global  flow  state  provides  a  challenging  problem.  Recent  research  in  fluidic  actuator 
development  (e,g,  synthetic  jets,  plasma  actuators,  vortex  generators)  provide  very  responsive  dynamics  such  that 
convective  time  delays  become  the  limiting  factor  in  terms  of  closed-loop  system  performance,  thus  model-free 
based  control  approaches  will  not  yield  the  best  performance. 

To  overcome  the  convective  time  delay  issue,  a  model  which  captures  the  dynamics  is  essential  for  predictive 
control.  Fluid  flow's  exhibiting  exploitable  instabilities  arc  governed  by  the  Navier-Stokes  equations,  a  set  of 
second-order,  non-lincar  partial  differential  equations.  Unfortunately,  the  complexity  of  the  governing  equations 
provides  little  to  no  insight  for  classical  or  modern  control  design  theory.  Also,  control  theory  is  very1  limited 
when  faced  with  high-dimensional,  extremely  non-linear  systems.  A  proven  wray  to  reduce  the  dimensionality  of 
the  problem  while  still  capturing  the  essential  system  dynamics  of  a  flow  problem  is  through  reduced  order  models 
(ROMs),  where  a  given  model  of  a  system  or  process  is  replaced  by  a  smaller  approximation  to  the  original  model. 
The  identified  ROM  can  then  be  used  as  a  predictive  model  to  minimize  the  time  delay  through  the  system. 

In  control  theory  time  delay  issues  are  a  common  problem  to  overcome.  Forms  of  internal  mode]  control  or 
model  predictive  control  can  be  used  to  accommodate  and  minimize  the  effect  of  such  issues.  In  particular,  Smith 
predictors  are  commonly  used  in  systems  with  pure  time  delays.  In  essence  the  Smith  predictor  is  broken  up  into 
two  loops,  one  of  which  computes  control  of  the  delay  free  model  using  a  primary  controller  (  typically  a  propor¬ 
tional  integral  derivative  control),  while  the  other,  minor  loop,  eliminates  the  actual  delayed  measurement;  thus, 
the  closed-loop  frequency  response  can  be  increased  beyond  the  performance  of  standard,  mod  el -free  approaches. 
The  downside  of  implementing  a  Smith  predictor  is  any  model  uncertainty  in  terms  on  non-linear  behavior,  inac¬ 
curate  time  delays,  or  varying  time  delays  drastically  hinders  the  performance  of  the  predictive  control. 

This  paper  uses  a  model  based  smith  predictor  to  control  a  non-linear  fluid  dynamic  problem.  The  predictive 
control  results  are  compared  to  that  of  the  open -loop  behavior  and  a  standard  model- free,  control  approach . 

A.  Flow  Field 

The  flow  field  around  a  slender,  axisymmetric  forebody  varies  dramatically  with  the  angle  of  attack.  Typically, 
four  flow  regimes  are  observed:  attached  flow  (0°  <  a  <  15°),  symmetric  vortex  flow  (15°  <  ft  <  40°), 
asymmetric  vortex  flow  (40°  <  ft  <  60°),  and  unsteady  wake- like  flow  {ft  >  60°)  [1  3],  Illustrations  of  these 
flow  regimes  arc  shown  Fig.  I ,  The  transition  between  the  symmetric  and  asymmetric  vortex  flow  is  due  to  an 
instability'  of  the  natural  flow  field.  Minor  geometric  imperfections  and  flow  perturbations  are  amplified  by  this 
convective  instability  and  divert  the  flow  field  away  from  the  symmetric  vortex  state  into  an  asymmetric  vortex 
state  in  wliich  either  port  or  starboard  vortex  will  separate  from  the  forebody  surface  and  cause  a  large  asymmetric 
pressure  distribution  on  the  ogive  surface. 


Once  this  occurs,  the  asymmetric  vortex  state  produces  a  significant  side  force  on  the  body,  typically  on  the 
order  of  40%  to  70%  of  the  normal  force  depending  on  angle  of  attack.  This  side  force  or  phantom  yaw  is  a 
detrimental  problem  for  the  stability  and  maneuverability  characteristics  of  slender  flight  vehicles  at  high  inci¬ 
dence  [4  6],  Numerous  experimental  and  numerical  investigations  for  a  range  of  geometric  shapes  and  Reynolds 
numbers  have  been  performed  to  understand  this  phenomenon. 

The  natural  tendency  of  the  flow  field  to  favor  an  asymmetric  vortex  stale  emphasizes  the  importance  of 
forebody  vortex  management.  A  variety  of  different  flow  control  methods  have  been  devised  for  different  forebody 
configurations.  Typically,  passive  control  techniques,  in  which  no  energy  is  added  to  the  flow,  make  use  of  nose 
strakes,  dimples,  nose  blunting,  and  other  geometric  modifications  to  reduce  the  flow's  instability  and  enforce  a 
symmetric  vortex  state  [4,7,8].  To  improve  vortex  management  performance  with  respect  to  passive  techniques 
and  open-loop  methods,  closed-loop  active  flow  control,  in  which  sensor  signals  are  used  to  estimate  the  flow 
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state  and  then  prescribe  the  actuation  input  to  reach  a  desired  state,  hold  the  promise  to  increase  stability  and 
maneuverability  characteristics  for  high  angles  of  attack  at  vary  ing  flow  conditions.  Active  closed-loop  flow 
control  presents  the  opportunity  to  not  only  regulate  side  force,  but  also  allow  for  the  possibility  to  attain  attitude 
control  of  a  slender  body  via  active  vortex  manipulation 

To  dale  very  little  research  has  been  performed  to  exploit  the  full  potential  of  feedback  flow  control  on  a\- 
isymmetric  forebodies,  Bernhardt  and  Williams  [9]  focused  on  regulating  the  side  force  over  a  range  of  angles 
of  attack  through  the  use  of  active  closed-loop  control  using  both  P1I>  and  artificial  neural  network-PD  control. 
The  overall  results  of  their  control  effectiveness  is  shown  in  Fig,  2a,  As  shown,  their  controller  was  able  to  effec¬ 
tively  manage  the  asymmetric  vortex  state,  all  but  eliminating  any  side  force  and  yawing  moment  up  to  an  angle 
of  attack  of  48°,  According  to  the  report,  the  highly  unstable  oscillatory  response  at  increasing  angles  of  attack 
exceeded  the  bandwidth  of  actuator  hardware  and  therefore  hindered  controller  perfo nuance  at  higher  angles  of 
attack.  Secondly,  Patel  et  al.  [ID]  were  also  able  to  successfully  close  the  loop  on  an  axisymmetric  forebody  by 
employing  deflectable  tabs  to  manipulate  the  vortex  state  and  achieve  a  desired  side  force  and  yawing  moment. 

1  lowcvcr,  the  use  of  deflectable  tabs  also  resulted  in  slow  actuation  changes  and  potentially  increased  the  drag  on 
the  model  (Fig.  2b). 


Figure  2,  Previous  experimental  closed -loop  tow  cool  ml  mu  Els  of  an  a  visym  metric  forefoody  at  high  angles  of  attack.  [9-12] 

Also,  closed-loop  experimental  results  from  the  US  Air  Force  Academy  [13]  which  implemented  plasma 
actuation,  surface  mounted  pressure  sensors  and  a  PID  control  algorithm  show  the  ability  to  reference  track  quasi¬ 
steady  side  force  set  points.  The  unsteady^  uctu  at  ions  (from  the  separated  shear  layer  along  the  bluff  body)  are 
unable  to  be  reduced,  as  shown  in  Fig  (??.  Because  plasma  actuators  have  very  responsive  dynamics,  the  lack  of 
the  ability  to  reduce  the  unsteady  fluctuations  proved  that  it's  not  associated  with  actuator  dynamics,  but  rather 
convective  time  delays,  as  the  closed-loop  bandwidth  was  only  limited  by  the  response  of  the  asymmetric  vortex 
state  due  to  an  introduced  disturbance  at  the  nose  of  the  ogive  and  propagation  time  of  the  disturbance. 
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Figure  X  Time  history  of  closed-loop  lest  with  the  reference  side  force*  Cyrrf%  in  dotted  black  and  estimated  side  force,  Cv>  shown  in 
hhtc,  |I3] 
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A.  Simulations 


II.  Setup 


The  geometry  considered  in  this  investigation  is  a  generic  von  Karman  ogive  which  is  a  part  of  the  Haack  series 
nose  cones  with  fineness  ratio  /r  =  3,5  and  a  diameter  of  D  =  0 ,  I  m,  The  nose  cone  shape  is  computed  from  Eq 
( 1)  and  (2),  where  the  origin  of  the  coordinate  system  is  at  the  nose  of  the  model. 


are  cos 


Hi  lit  2fl„) 


A  cylindrical  aft  body  with  length  Laff  =  0.05  m  is  added  lo  the  end  of  the  ogive  forebody.  Figure  4  shows  the 
geometry  used  in  both  the  experiment  and  simulation. 

The  CFD  solver  used  is  Cobalt,  an  unstructured  finite- volume  code  developed  for  the  solution  of  the  compress¬ 
ible  Navier-Stokcs  equations.  The  basic  algorithm  is  described  in  Strang  et  al.,  although  substantial  improvements 
have  been  made  since  then  [14].  The  numerical  method  is  a  cell-centered  finite  volume  approach  applicable  to 
arbitrary  cell  topologies  (c.g,  hexahedra*  prisms*  tetrahedra).  The  spatial  operator  uses  a  Riemann  solver*  least 
squares  gradient  calculations  using  QR  factorization  to  provide  second  order  accuracy  in  space.  A  point  implicit 
method  using  analytic  first-order  inviscid  and  viscous  Jaeobians  is  used  for  advancement  of  the  discretized  sys¬ 
tem.  For  time-accurate  computations*  a  Newton  sub-iteration  scheme  is  employed,  resulting  in  a  method  that 
is  formally  second  order  accurate  in  time.  For  parallel  performance,  Cobalt  utilizes  the  domain  decomposition 
library  ParMETIS  to  provide  optimal  load  balancing  with  a  minimal  interface  between  zones  [15]. 

To  compute  the  flow  around  the  model,  an  unstructured  grid  was  generated  which  contained  16  million  total 
elements.  Near  the  surface  of  the  body,  a  boundary  layer  grid  wras  used  to  provide  a  dense  spacing  of  grid  points 
near  the  model  due  to  the  sensitivity  of  the  flowfield  to  small  disturbances.  A  spherical  farfield  boundary  was 
placed  40  model  diameters  away  from  the  origin  to  minimize  the  influence  of  possible  pressure  reflections.  A 
schematic  of  the  model  geometry,  used  in  the  accompanying  experiment,  and  replicated  computationally  is  shown 
in  Fig.  4. 


liaised  RScctrodc 
<jf  Phiisintp  Acluator 


Pressure  Transducer 
x/D  ^2.0  m  ±ao 


Pressure  Transducer 

x/d  =  3,0  =  ±m( 


Figure  4.  Schematic  cl  riming  of  experimental  axisyinmetric  forebotly.  The  pressure  transducer  locations  ]I6|  and  plasma  plasma 
actuators  arc  shown. 


The  experiment  employed  two  single  dielectric  barrier  discharge  plasma  actuators  at  the  tip  of  the  ogive  lo 
augment  the  flow  field  as  described  in  Faglev  et  al  [?]♦  To  replicate  this  actuation  mechanism  within  the  CFD 
simulation,  two  user- specified,  mass-blowing  patches  vvere  placed  at  ±90°  from  the  windward  meridian,  1  cm 
downstream  from  the  lip  of  the  model.  These  patches  were  2  cm  long  with  a  width  of  0.015  cm.  The  mass- 
blowing  patches  were  set  to  blow  at  6  =  30IJ  (Fig.  5a)  to  mimic  the  plasma  actuators  in  the  accompanying 
experiment  [17].  Finally,  a  small  pin  shaped  disturbance  was  added  at  90°  to  the  starboard  side  of  the  model. 
The  pin  had  a  diameter  of  0,01  cm  and  was  0.0045  cm  tall.  The  center  of  the  pin  was  placed  0.4  cm  downstream 
from  the  tip  of  the  model.  The  pin  provided  an  initial  geometric  disturbance  to  initiate  a  natural,  deterministic 
asymmetric  vortex  state.  Figure  5a  shows  a  zoomed  in  view  of  the  tip  of  the  model,  showing  the  pin  and  the 
actuator  design. 

Surface  mounted  pressure  transducers  were  placed  in  the  ogive  to  estimate  the  instantaneous  side  force,  be¬ 
cause  realtime  force  measurements  are  susceptible  to  frequency  responses  which  are  limited  to  the  inertia  of  the 
body  and  are  not  feasible  in  actual  application.  The  locations  of  the  pressure  transducers  are  shown  in  Figure  4. 
The  sensor  placement  was  determined  through  an  optimal  genetic  algorithm  as  documented  by  fagley  et  al.  f  i  6]. 
This  arrangement  of  sensors  minimized  the  linear  estimation  between  the  four  pressure  transducer  measurements 
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and  side  force  measurements  in  forced  and  unforced  CFD  simulations.  The  estimated  force  is  computed  by  the 
linear  model, 

Cy(t)^CP(xs,t).  (3) 

where  P(xlsJ)  is  the  array  of  pressure  measurements  at  surface  locations  x C  is  the  vector  of  coefficients 
which  map  pressure  measurements  into  the  estimated  side  force,  Cy.  The  port  side  plasma  actuator  is  also  shown 
in  Figure  4.  Detailed  information  on  the  fabrication,  implementation  and  open  loop  performance  of  the  SDBD 
plasma  actuator  can  be  found  in  [18]. 


u)  Stiirlmard  and  port  mass  blowing  slots  b)  Cross-seetionaJ  view  of  the  grid  at  X/D  =  2. 


Figure  5.  Starboard  and  port  mass  blowing  slots  shown  in  (a)  along  with  front  (b)  view  of  the  computational  grid  used.  The  grid 
contains  6, 711, 157  nodes  and  15, 568, 273  dements. 


111.  Results 


All  of  the  figures  and  data  presented  are  normalized  by  the  convective  flow  through  time  for  the  free  stream 
velocity  as  given  by, 

tU°°  (4) 


with  the  corresponding  convective  frequency. 


fr  = 


L  ' 
L 


(5) 


tux 

To  characterize  the  effectiveness  of  the  actuator  within  the  CFD  simulation  the  non-dimensional  momentum  coef¬ 
ficient  given  by. 


„  \JApV{V.h)dA\ 

PooUlA  * 


(6) 


was  used. 


A,  Approach 

To  reduce  the  effect  of  convective  time  delays,  the  closed-loop  response  time  is  minimized.  For  this  a  modified 
smith  predictor  is  implemented  to  increase  the  crossover  frequency  and  phase  margins  of  the  closed-loop  system. 
The  chosen  control  structure  is  shown  in  Figure  6.  This  controller  uses  □  linear  model  without  the  pure  time  delay 
as  a  future  predictor  of  the  plant  response  to  improve  upon  controller  responsiveness.  The  measured  output  is  then 
negated  by  the  delayed  model  response,  thus  if  the  plant  is  exactly  equal  to  the  model  the  pure  time  delay  can  be 
accounted  for.  The  control  algorithm  or  compensator  (G,  (s)),  is  typically  a  proportional  -  integral  -  derivative 
algorithm.  The  advantage  is  such  that  the  error  term,  or  control  input,  now  operates  on  an  estimated  model 
response  at  a  convective  flow  through  time  advancement. 

To  implement  this  type  of  predictive  control,  a  model  is  now  necessary  of  the  underlying  physical  dynamics. 
The  system  response  was  modeled  using  a  linear  system  parameterization.  The  input-output  relationship  for  this 
system  is  given  by, 

Y(s)  =  G'ta)U(s).  (7) 

The  structure  of  the  model  in  continuous  time  has  the  form. 


G  (  \  =  —  K  +  am * 1  1  +  Q™—2*m  2  ■  -  a  t  -H  +  at j 

*  *  D(&)  6  sn  +  \Su  l  H-  K-2S71  ~2  ■  ■  &!  s  +  &0 


(8) 


for  a  linear  system  with  m  zeros,  n  poles,  and  a  pure  time  delay,  e&5.  A  pure  time  delay  is  necessary  because  of 
the  physical  convective  time  from  which  a  disturbance  is  introduced  to  the  flow  and  reaches  the  sensing  location. 


5  of  1A 


American  Institute  of  Aeronautics  and  Astronautics 


d 


Figure  6*  Smith  predictor  control  structure  for  control  of  asymmetric  state  on  a  von  Karman  Ogive 


Different  system  identification  techniques  exist  for  parameterizing  suitable  models  G(s),  determining  the  poly¬ 
nomial  orders  m  and  n,  and  solving  for  the  coefficients  of  the  polynomials  in  the  numerator  and  denominator. 
The  technique  for  lime  domain  identification  used  in  this  effort  is  the  Prediction  Error  Method  (PEM),  The  PEM 
technique  is  an  iterative  approach  to  determining  a  minimum  error,  in  a  least  squares  sense,  for  a  selected  model 
structure  given  by  an  autoregressive  moving  average  (arma)  system.  This  approach  works  m  the  continuous  and 
discrete  time  domains  and  is  applicable  for  multi-input  multi-output  time  domain  data  sets.  The  PEM  technique 
directly  identifies  with  a  model  structure  of  the  form  of  Rq.  (8).  In  particular,  a  model  structure  as  shown  in  (9) 
was  suitable  for  capturing  the  temporal  dynamics  due  to  a  step  input. 


Gs(s)  =  Kp 


i  4- 

(1  4-  2{Tws  4  (T^s)2) 


e 


-Tds 


(9) 


Because  high-fidelity  CFD  simulations  are  CPU  demanding  and  take  on  the  order  of  60  hours  fora  single  sim¬ 
ulation,  this  does  not  allow  for  an  amenable  environment  for  control  algorithm  design  and  tweaking.  Therefore, 
an  intermediate,  non-linear  model  must  be  developed  to  simulate  the  closed- loop  dynamics.  Previous  results  show 
that  three  forcing  regimes  exist.  At  lower  forcing  magnitudes  of  the  jet  actually  create  a  separation  bubble 
which  accelerates  the  flow,  delaying  separation  and  causing  an  attached  vortex.  At  medium  forcing  magnitudes, 
the  jet  induces  separation  causing  a  detached  vortex,  and  at  highest  forcing  magnitudes,  the  jet  produces  a  coanda 
like  effect  and  causes  an  attached  vortex  state  [19].  Thus,  for  the  purposes  of  this  paper,  only  forcing  magnitudes 
which  induce  separation,  or  induce  a  detached  vortex  arc  considered  because  the  largest  range  of  side  force  mag¬ 
nitudes  are  attained  in  this  forcing  regime.  In  this  forcing  regime,  although  the  temporal  dynamics  due  to  a  step 
input  are  relatively  constant,  the  steady-state  response  of  side  force  is  non-linear  as  shown  below  in  Fig  8a, 

A  Ham  me  rstein- Wiener  model  structure  {NLHW)  was  selected.  This  model  structure  applied  a  nonlinear 
function  to  both  the  inputs  and  outputs  of  a  linear  model  The  NLHW  model  was  built  using  sigmoidnets  to 
account  for  the  input  nonlinearity  and  wavenets  to  account  for  the  output  nonlinearity.  Furthermore,  the  model 
was  trained  using  only  the  open-loop  data,  and  validated  against  the  closed-loop  data.  To  determine  the  number  of 
sigmoidnets  and  wavenets  needed,  the  model  was  trained  allowing  each  to  vary  from  one  to  ten,  and  the  average 
fit-error  across  all  the  data  sets  was  calculated,  including  the  closed-loop  data.  The  final  results  indicated  that 
eight  sigmoidnets  and  lour  wavenets  provided  the  best  fit  to  the  data.  This  creates  a  black-box  model  structure 
built  solely  for  the  purposes  of  data  replication,  such  that  it  provided  an  accurate  model  for  controller  development 
The  results  of  a  step- input  in  forcing  as  well  as  the  steady-state  response  are  also  shown  in  Fig.  ??, 


B,  Transient  Dynamics 

To  model  the  system  dynamics  of  the  asymmetric  vortex  state  behind  the  von  Karman  ogive  at  high  angles  of 
attack,  open- loop  CFD  simulations  were  conducted  to  understand  system  character i sties  in  tonus  of  stability/bi¬ 
stability,  controllability,  observability,  and  linear/ non-linear  behavior.  The  flow  behind  an  axisymmetric  slender 
body  has  previously  been  shown  to  be  completely  bistable  at  a  sufficiently  large  angle  of  attack  and  small  Reynolds 
numbers.  Because  few'  tests  have  been  conducted  at  a  Reynolds  number  in  the  range  of  the  current  simulation 
(Re  =  156,000)  with  this  unique  geometry,  determining  if  a  bi-stable  or  proportional  flow  regitne  exists  is  critical 
in  designing  a  suitable  model  structure  and  control  system  design. 

The  open-loop  dynamics  are  determined  purely  from  step  and  impulse  actuation  inputs.  The  step  inputs 
were  varied  at  different  actuation  inputs  in  terms  of  plasma  voltages  and  blowing  coetFie  rents  for  the  experiments 
and  simulations,  respectively.  The  impulse  input  was  varied  over  a  range  of  flow  through  limes.  Figure  7a 
shows  the  comparison  of  the  step  response  for  the  experimental  measurements  of  the  side  force,  the  numerical 
simulation  and  the  model  estimation.  Both  of  the  numerical  and  experimental  data  were  offset  so  that  the  initial 
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state  was  zero.  As  shown,  the  agreement  is  very  good  among  experimental,  numerical  and  model  response.  The 
step  input  data  ensemble  served  for  the  training  data  for  which  the  linear  model  was  fit.  Figure  7b  shows  the 
experimental,  numerical  and  model  responses  to  an  impulse  response  which  varied  with  duty  cycle  by  percentages 
of  the  normalize  flow  through  lime.  The  impulse  dam  set  served  as  the  validation  for  the  model  verification.  Very 
good  agreement  between  all  three  systems  existed. 


Figure  7.  a)  Step  response  for  experiment  |btue|,  simulation  j  green],  and  LTI  model  prediction  |modelf  b>  Impulse  response  for 
experimeutl-J,  simulation!-!,  and  LTI  model  M  at  varying  percentages  of  r. 
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a)  Open-loop  summary 
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Figure  8.  a)  Open-loop  su  mmary  Tor  steady  state  response  of  Cy  due  to  step  input  of  Op  (positive  starboard  forcing  and  negative  port 
forcing),  Linear  and!  non-linear  representations  also  plotted,  b)  Conformal  mapping  of  Op  — *  C*t  using  a  nil  ion  ill  polynomial  best  fit. 


C.  Closed-Loop  Results 

The  LTI  model  developed  above  provides  a  relationship  between  the  input  momentum  coefficient  and  the  resultant 
estimated  side  force.  To  compare  the  results  of  a  predictive  control,  a  simple,  model-free  controller  was  initially 
developed.  A  proportional-integral  controller,  G?c{a)  —  Kp  +  was  selected  to  close  the  loop. 

such  that  the  dynamics  of  a  closed- loop  controller  can  be  evaluated.  A  schematic  of  the  feedback  system 
implemented  is  shown  in  Fig.  ??  and  formulated  in  Eq.  (10), 

Y=  [  l  +  c2\  l  +  C.C,  ] 

where  r  and  d  are  the  reference  and  disturbance  inputs,  respectively. 

With  the  system  shown  in  Eq.  ( 10} ,  the  transfer  function  between  different  input-output  pairs  can  be  analyzed 
for  various  forms  of  GTr(.s).  As  an  initial  control  design,  a  standard  PI  controller  for  <7,(  *0  was  implemented, 

Ge(«)  =  Kp+  ^  (If) 

where  Kp  and  A,  arc  the  proportional  and  integral  terms,  respectively.  No  derivative  gain  was  used  since  the 
derivative  term  resulted  in  an  instability  in  the  transfer  function,  GSGCJ{  1  4-  G&Gr),  purely  due  to  the  time  delay 
in  the  system.  The  gains,  A\  and  A'p,  were  selected  to  achieve  high  fluidic  response  times  while  maintaining 
adequate  stability  margins.  As  shown  in  Fig.1??,  the  open-loop  phase  margin  is  05°  and  the  open-loop  gain 
margin  is  1 .63  or  5  dB. 
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Figure  9. 


IV,  Conclusion 

The  asymmetric  vortex  regime  of  a  axisymmetric  bluff  body  with  a  fineness  ratio  of  3.5  is  experimentally 
and  computationally  studied  at  a  Reynolds  number  of  1 56,000.  Both  port  and  starboard  plasma  actuators  are 
used  to  introduce  fluidic  disturbances  at  the  lip  of  the  ogive  in  the  experimental  setting.  Mass  blowing  ports 
which  introduced  momentum  in  the  flow  were  used  in  the  simulations*  These  disturbances  are  amplified  through 
the  flow's  convective  instability  to  produce  a  deterministic  port  or  starboard  asymmetric  vortex  state  {i.e.  side 
force).  Accurate  control  or  manipulation  of  this  asymmetric  vortex  phenomenon  holds  the  potential  for  increased 
maneuverability  and  stability  characteristics  of  slender  flight  vehicles* 

Unforced  and  open-loop  experimental  tests  were  carried  out  to  understand  and  quantify  the  vortex  dynamics. 
Step  and  impulse  inputs  provided  the  necessary  dynamics  and  diverse  training  and  validation  data  sets  for  the 
formulation  of  a  linear  time  invariant  dynamical  model.  Standard  linear  system  identification  approaches  were 
implemented  to  represent  the  training  data  set.  In  particular  prediction  error  was  used  to  capture  the  asymmetric 
vortex  dynamics.  These  methods  were  validated  by  time  and  frequency  domain  methods.  The  measurements  and 
modeling  methods  showed  the  cutoff'  frequency  of  the  flow  to  be  around  two  flow  through  times,  i.e*  the  time  it 
takes  a  particle  to  flow  from  the  tip  of  the  model  to  the  base  of  the  ogive  section. 

A  closed- loop  system  was  designed  such  that  the  unforced  fluid  dynamics  and  measurement  noise  were  mod¬ 
eled  as  an  output  disturbance.  The  prediction  error  model  was  well  suited  for  this  system.  A  PID  controller  was 
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Figure  II. 


Non-dimensional  time,  r 

Figure  12. 


implemented  in  the  closed  loop  system  and  designed  for  adequate  disturbance  rejection  and  reference  tracking  per¬ 
formance.  The  closed  loop  transfer  functions  were  analyzed.  A  time  simulation  was  shown  in  which  the  controller 
was  able  to  guide  the  asymmetric  vortex  state  to  an  arbitrary  asymmetric  pressure  distribution  while  adequately 
regulating  the  disturbances.  The  control  was  then  validated  in  both  experiment  and  Navier  Stokes  numerical  sim¬ 
ulations.  Very  good  agreement  existed  between  the  two  systems.  The  dynamics  showed  that  frequencies  larger 
than  a  flow  though  lime  could  not  be  controlled.  This  is  mainly  due  to  the  actuator  placement  at  the  nose  of  the 
axisymmetric  body  and  the  large  convective  time  for  a  disturbance  to  reach  the  end  of  the  geometric  body.  Also 
it  is  shown  that  the  high  frequency  disturbance  is  mainly  due  to  the  separated  shear  layer  along  the  side  of  the 
axisymmetric  body.  The  modeling  and  understanding  of  the  interaction  between  these  primary  vortices  and  the 
separated  shear  layer  is  critical  for  increasing  the  performance  of  flow  control  of  an  ax  isy  metric  body  at  high  angle 
of  attack. 
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